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The dynamics of gluons and quarks in a relativistic nu- 
clear collision are described, within the framework of a classi- 
cal mean-field transport theory, by the coupled equations for 
the Yang-Mills field and a collection of colored point parti- 
cles. The particles are used to represent color source effects 
of the valence quarks in the colliding nuclei. The possibil- 
ities of this approach are studied to describe the real time 
evolution of small x modes in the classical effective theory in 
a non-perturbative coherent manner. The time evolution of 
the color fields is explored in a numerical simulation of the 
collision of two Lorentz-boosted clouds of color charged par- 
ticles on a long three-dimensional gauge lattice. We report 
results on soft gluon scattering and coherent gluon radiation 
obtained in SU(2) gauge symmetry. 



I. INTRODUCTION 

Experiments with relativistic heavy ions at center-of- 
mass energies reaching 100 GeV/u will soon search for 
a new phase of nuclear matter [|l]j2|. It was argued that 
the extraordinarily high energy and particle number den- 
sities reached in central nuclear collisions at RHIC ||] 
could lead to rapid local thermalization of matter [Q and 
thus to the formation of the so-called quark gluon plasma 

i- 

One of the theoretical challenges in this context is the 
development of a description on the basis of quantum 
chromodynamics (QCD) of the processes that may lead 
to the formation of the locally equilibrated quark gluon 
plasma state in these nuclear reactions. Predictions on 
whether and how this state could be formed in a nuclear 
collision depend in a sensitive way on the initial con- 
ditions. At present, our understanding of these condi- 
tions is rather poor. Perturbative parton-parton scatter- 
ing precedes local thermo-dynamical equilibrium in the 
early stage and may partly determine the dynamics of 
thermalization. Indeed, event generators based on indi- 
vidual parton-parton processes Q have been developed. 
However, their predictions differ widely. To improve our 
understanding of the early stage, we have to consider the 
fact that the initial state is described by coherent parton 
wave functions and glue fields. Consequently, coherent 
multiple scattering must play a role in the early stage of a 
collision. The simplest approach to this picture would be 
to start with a mean field description of the color fields 



in an effective theory while hard modes are described 
through classical particles representing partons. 
The possibility of a description of inelastic gluon pro- 
cesses by means of the nonlinear interactions of classical 
color fields has been proposed some time ago . Some 
years ago, this scenario was examined in studies of the 
collision of two transverse polarized Yang-Mills field wave 
packets on a one-dimensional gauge lattice Q . These cal- 
culations showed that the interaction between localized 
classical gauge fields can lead to the excitation of long 
wavelength modes in a way that is reminiscent of the 
formation of a dense gluon plasma. Recently, a similar 
study has been carried out for collisions of Yang-Mills 
wave packets on a three-dimensional gauge lattice |10|. 
Here, we go a step further and add color charged particles 
to the system acting as sources of the color fields. Instead 
of two wave packets, we collide two clouds of such par- 
ticles which shall describe the color charge distributions 
generated by the valence quarks in two colliding nuclei. 
The wave packets are now replaced by the colliding gluon 
fields which are generated by classical color sources while 
they propagate through the lattice, representing the fast 
moving nuclear valence quarks. The inclusion of particles 
also permits a comparison with perturbative calculations 
of gluon radiation from colliding quarks 1 1 1 13] . 



We address the question whether indications for the for- 
mation of a quark gluon plasma can also be observed in 
such simulations. For simplicity, we leave out the hard 
collisions between the particles and focus first only - ex- 
cept from the time dependence of the color sources - on 
the soft dynamics in the gluon fields which are described 
on the lattice and correspond to modes with a small light 
cone momentum fraction x. We further leave out the 
hard modes of the gluon fields which may also be rep- 
resented by particles. The particles may in principle be 
used in the future to include the semi-hard dynamics at 
larg er x following the idea of the parton cascade model 
||l3| . While this "probabilistic" approach provides a rea- 
sonable description of large transverse momentum pro- 
cesses at large x, QCD coherence effects become impor- 
tant as we go to small x or alternatively, towards central 
rapidities [Q. What is needed therefore to describe the 
collision of the "wee" partons is a wave picture where 
coherent multiple scattering is fully taken into account. 
At the limit of extremely high collision energies, the 
model presented here has in principle to converge against 
a random light-cone source model which was initially pro- 
posed by McLerran and Venugopalan |jl^ and later fur- 
ther developed by these authors in collaboration with 
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Ayala, Jalilian-Marian, Kovner, Leonidov, and Weigert 
|pi|-[19|. This model may therefore be used as a refer- 
ence in the hmit u/c — > 1. The McLerran-Venugopalan 
model considers very large nuclei moving with nearly the 
speed of light which consequently appear in the labora- 
tory frame as "pancakes" of almost zero thickness in the 
transverse plane. The classical gluon field at small x for 
a nucleus in the infinite momentum frame is obtained 
by solving the Yang-Mills equations in the presence of a 
static source of color charge on the light cone. In our ap- 
proach, the solution of the equations of motion requires 
one to take the longitudinal extension of the nucleus into 
account, i.e. one must not take the nucleus to be in- 
finitely thin in the longitudinal direction as assumed in 
1^ . Following Kovchegov and Rischke |^ , we therefore 
take the finite longitudinal extension of the nuclei into 
account. This of course requires that the velocities of the 
nuclei are slightly less than the speed of light. 
A major goal of the description, presented in this paper, 
is to keep the possibility open for an inclusion of parton 
cascades in the further development while the coherent 
physics of the random light-cone source model at small x 
is described at the same time. This goal can only be ac- 
complished in several steps. The approach presented here 
contains no restriction to central rapidities and provides 
therefore an extension to larger x taking into account 
semi- hard gluons. 

In section 2, we present a brief description of the model 
and the formulation on a SU(2) gauge lattice. In section 
3, we discuss the preparation of the initial state of the 
nuclei and show numerical results. In section 4, we evolve 
the collision and present results for the time evolution of 
the glue fields in the collision. A summary and conclusion 
is given in the section 5. 



II. DESCRIPTION OF THE MODEL 

A set of equations describing the evolution of the phase 
space distribution of quarks and gluons in the presence 
of a mean color field, but in the absence of collisions, 
was proposed by Heinz pl| , p^ . This non-Abelian gen- 
eralization of the Vlasov equation can be considered as 
the continuum version of the dynamics of an ensemble of 
classical point particles endowed with color charge and 
interacting with a mean color field. Denoting the space- 
time positions, momenta, and color charges of the parti- 
cles by , and g", respectively, where i = 1, . . . , iV is 
the particle index, the equations p3| for this dynamical 
system read: 



The moving particles generate a color current J'^ = 
J^T°/2 which forms the source term of the inhomoge- 
neous Yang-Mills equations 

[V^,T^''{x)] =gj''ix) 

= 9y^Q^{tf^5{x^xm, (2.4) 

i 

for the mean color field. These equations have recently 
been used in the weak-coupling limit {g -^V) to simulate 
the effects of hard thermal loops on the dynamics of 
soft modes of a non-Abelian gauge field at finite temper- 
ature [||,|6|. 

Here wc use Eq. (pj|-2.4) to describe the interactions 
among the soft glue field components of two colliding 
heavy nuclei. Transverse modes are included by using a 
3-dimensional spatial gauge lattice. The short-distance 
lattice cut-off a separates the regime in transverse mo- 
mentum where the dynamics of gluons is perturbative 
(large /ct) from that where perturbation theory fails 
(small fcx) and which is of interest to us here. The inter- 
action with the mean color field allows for an exchange 
of an arbitrary number of soft gluons. In combination 
with a parton cascade, the screening of the soft compo- 
nents of the gauge field by perturbative partons 28 
is taken into accou nt natur ally by the nonlinear nature 
of the coupled Eq. (2.1-2^). This will be discussed in a 
separate publication |29|] . 

We represent the valence quarks of the two colliding nu- 
clei as point particles moving in the space-time contin- 
uum, and interacting with a classical gauge field defined 
on a spatial lattice in continuous time. Here we ne- 
glect the collision integrals describing hard interactions 
between the particles. In the spirit of the statistical na- 
ture of the transport theory, we split each quark into a 
number Uq of test particles, each of which carries the frac- 
tion qo = Qo/ riq of the quark color charge Qo = ^^3/4. 
For the gauge group SU(2) adopted here, each nucleon 
is represented by two quarks, initially carrying opposite 
color charge. A lattice version of the continuum equa- 
tions (|2.1|-]2.4D is constructed ^ by expressing the field 
amplitudes as elements of the corresponding Lie algebra, 
i.e. A^,J-^,y,£k,Bk G LSU(2). As in the Kogut-Susskind 
model of lattice gauge theory |^ we choose the temporal 
gauge ^0=0 and define the following variables. 



Ux.i = exp(-i5ai^i(x)) 

UxM = Ux,k Ux+k,l ^x+k+l,-k l^x+l,-l 

Consequently, we have 
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for the electric and magnetic fields, respectively. The 
lattice constant in the spatial direction j is denoted by 
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ttj. As one can see from (2.5), the gauge field is expressed 



in terms of the link variables U^.i e SU(2), which represent 
the parallel transport of a field amplitude from a site x 
to a neighboring site {x + l) in the direction I. We choose 
l^x,k and £x,k as the basic dynamic field variables and 
solve the following equations of motion for the fields. 



Ux,k{t) = igak £x.k{t)^x,k{t) 



(2.8) 



£x,k{t) — 
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- ^i-uW {i^x-i,kiit)-i^Li,ki(t)) ^-mW } 

- gJxAt)- (2-9) 

In the time evolution of the particle lattice system, these 
equations have to b e sol ved simultaneously together with 
the equations (2J) to ( ^.3[ ). The dynamics of the par- 
ticles is described within a dual lattice which is defined 
by the edges of the Wigner-Seitz cells of the original lat- 
tice. The charge of all particles in the volume of one cell 
is associated with the center of the cell represented by 
a point of the real lattice. As long as a particle prop- 
agates in the cell its charge vector is not changed. If a 
particle crosses the boundary of a dual box, the charge is 
transported parallel to the center of the neighbor cell into 
which the particle enters. Exchange of energy between a 
particle and the field occurs only at the transition from 
one dual box into another. The particle is reflected from 
the wall if its kinetic energy is smaller or equal than the 
work required to move the charge from the old box into 
the new box. In the other transition into the 

new box occurs where the amount of energy required to 
make the transition is either subtracted or added to the 
kinetic energy of the particle depending on the sign of 
the product Tr(Qi£x,i). This mechanism provides a cou- 
pling between the fields and the particles which allows to 
generate color fields through color charge currents. The 
method has the advantage that it exactly preserves the 
law of Gauss. It also preserves the conservation of energy 
which is thus only violated through numerical errors of 
the order At^ in the time integration of the field equa- 
tions (2.8) and (2.9). On the other hand, the method. 



as simple as it is, is connected with serious limitations 
as for example the violation of causality when the time 
step width At is smaller than the lattice constants o;. 
At can not be chosen larger than o; b ecau se the time 
integration of the equations (2.8) and (2.9) is instable 
for At > min(oi : I = 1,2,3) These limitations do not 
yet allow a study of the dynamics of the particles in a 
microscopically correct way. Artifacts, as for example a 
trapping of particles of too small kinetic energy in dual 
cells can occur. For propagating nuclei, the field energy 
associated with a single link can occasionally exceed the 
kinetic energy of a particle resulting in the emission of 
particles due to backwards scattering from the link. This 
corresponds to the (unphysical) emission of nuclcons from 
the moving nuclei, which we want to avoid. 



However, a slight modification of the above described 
procedure for the wall transition of particles allows to use 
the particles to simulate color charge fluctuations and to 
generate coherent color fields: We leave the kinetic en- 
ergy of a particle unchanged whenever it crosses a cell 
boundary of the dual lattice and only transport its color 
charge to the neighboring lattice cell. This is equivalent 
t o se tting the Coulomb force in the r.h.s. term of Eq. 
( ^.3[ ) to zero. In this way, the collection of particles acts 
as a classical color source propagating through the lattice 
with constant velocity. 



III. THE INITIAL STATE 

In the following, we briefiy outline how, in principle, a 
simulation of a nuclear collision is performed within the 
model described in the previous section. Mainly, how- 
ever, we discuss the question how to generate the "initial 
state" . 

For a central collision of two Pb nuclei one would expect 
that 34 nucleons in a row collide on the collision axis (z- 
axis) . This number is desired from the following estimate. 
We denote the nucleon radius by r„, the proton radius 
of a nucleus by i?p and the neutron radius of a nucleus 
by i?7v The average number of nucleons in a cylinder of 
the radius r„ around the collision axis is determined by 



R 



N 



N 



(3.1) 



1.32 fm, Rp ^ 5.47 fm, Rn = 
we find the result n = 17.11 for 



With the values 
5.80 fm from Ref. i 

In accordance with the spatial extension of a nucleon 
we choose a lattice extension of 1.2 fm into both trans- 
verse directions. The lattice spacing is taken o/ = 0.3 
fm in each direction, thus 4^ x 40 lattice points cover 
the volume of 17 nucleons in one row. More correctly, 
since there is no Pauli-Principle between neutrons and 
protons there are two superposed rows of nucleons. One 
row contains 8 protons and the other 9 neutrons. These 
details, however, are far beyond the reach of our cal- 
culation, and we simply take 16 nuclcons or 32 quarks 
(in SU(2)) in our initial state. We also use the rounded 
value r„ — 1.2 fm. The coverage of two complete nu- 
clei in the transverse plane requires much larger lattices 
and remains a challenge for more extended calculations 
in the future. The lattice is closed to a 3-torus which 
means that we apply periodic boundary conditions in all 
three spacial directions. The fact that we use a lattice 
with transverse extensions much smaller than the radius 
of a nucleus leads to restrictions. The first of these is 
that we can not simulate the full transverse dynamics 
of a collision. Such simulations for central collisions are 
possible but they require lattices with transverse exten- 
sions at least twice as large as the radius of a nucleus. 
Second, for lattices with small transverse extensions, the 
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periodic boundary conditions in the transverse directions 
define a situation which is equivalent to the assumption 
that the nuclei have infinite extension into transverse di- 
rections. Therefore, in this respect, our description stays 
firmly within the framework of the model of McLerran 
and Vcnugopalan Larger lattice sizes will be nec- 

essary to carry out a meaningful Fourier analysis of the 
fields in transverse directions and to compare results of 
both models. However, the transverse extension of a nu- 
cleon can be covered at least to describe the initial state 
of nucleons. 

As already mentioned above, a dual lattice is superim- 
posed on the original lattice in such a way that the lattice 
points are located in the centers of the cells of the dual 
lattice. The dual lattice cells are used to associate parti- 
cles with lattice sites and thus to define the source terms 



in (2^). To generate the initial configurations of the nu- 
clei, we randomly distribute color charged massless parti- 
cles over the volume such that each lattice cell is occupied 
by an even number Ub of particles. The total initial color 
charge is zero in each box corresponding to a neutral 
color charge distribution. Momenta with opposite but 
random directions and Boltzmann distributed absolute 
values are assigned to each pair of particles Q The Mo- 
menta are renormalized such that the total initial kinetic 
energy of the particles in one nucleus agrees with half of 
the correct mass of 16 nucleons {M = 8 • m„ = 8 • 939 
MeV) while we assume that half of the nuclcon energy is 
carried by glue fields. This procedure fixes the tempera- 
ture parameter T in the Boltzmann distribution for the 
particles. 

The initial configuration of a single nucleus in its rest 
frame is o b tain ed through t he evolution of the equa- 
tions {2A ~ 2^), (2^), and ( ^ ) over a long period of 
time starting from the initial fields £x,k{t = 0) = 0, 
Ux,k{t = 0) = I2. The simulation of a colhsion requires a 
Lorentz boost of each nucleus into the center of velocity 
frame of both nuclei. In the example presented below, 
the kinetic energy is 100 GeV/u for which 7 = 106.5. 
Both nuclei are mapped into their initial positions on a 
large Lorentz-contracted lattice right after the boost of 
particle coordinates Xi,pi and field amplitudes. 
At this point, it has to be emphasized that it is not pos- 
sible to boost fields on a lattice. We have extensively 
explored the possibility to generate the color fields be- 
fore the Lorentz-boost. In this case one would generate 
the initial state of the nucleus before the boost by evolv- 



^ In the "real" initial state of a nucleus, each nucleon of any 
shell is in its ground state with a temperature T = Q. On the 
other hand each nucleon has a finite energy according to its 
rest mass. Only a full quantum theoretical description can 
lead to a finite ground state energy at zero temperature. In 
our classical approach T plays the role of a parameter in the 
momentum distribution employed, to adjust the ground state 
energy. 



ing the equations of motion in time until an equilibrium 
between kinetic particle energy and field energy is ac- 
complished. This would require to employ the algorithm 
without the modification mentioned in section 2. The 
Lorentz-transformation of such a "initial state" , however, 
does not lead to a configuration which propagates stably 
along the lattice. One reason is that the lattice disper- 
sion relation is not boost invariant. After the boost, the 
Fourier components of the fields do not match with the 
finite discrete Fourier spectrum of the lattice. Another 
reason is that the Hamiltonian of the system does not 
commute with the boost operator. We have found it 
most practical to boost the particle coordinates a;f,pf 
at time t — Q for both nuclei, wherafter the fields are 
generated from the initial conditions £x.k(t = 0) = and 
Ux,k{t = 0) = I2. It has been argued [Q that the coupling 
constant as should be taken at an effective momentum 
scale on the order of 2ttT corresponding to 5 = 2. Start- 
ing from 5 = 0, the coupling constant is increased slowly 
until it reaches a final value g — 2. This method leads 
to quasi-stable configurations propagating over distances 
which extend over several thousand lattice points in the 
longitudinal direction. It avoids the problem with the 
lattice dispersion since the discretized k-space is fixed 
after the boost. The adiabatic increase of the coupling 
constant avoids the excitation of instable high frequent 
modes in the color fields . 

However, it requires large initial distances between the 
nuclei because the charge fluctuation converges slowly 
against a final constant value before the collision. Also, 
the fields should carry enough energy. For nucleons with 
the correct mass (m„ = 939 MeV) about half of the total 
energy are carried by glue fields. If we want to satisfy this 
condition to a better approximation then also semi-hard 
gluons have to be taken into account. The largest fraction 
of these 50% is carried by soft gluons with a Bjorken scale 
X < 0.1, followed by semi-hard gluons with 0.1 < x < 
0.2. In the original work of McLerran and Vcnugopalan 
|l5| , only the valence quarks were taken to be sources 
of color charge which gave /i^ ~ A^^^fm~^. Hence fi ^ 
^QCD only for nuclei much larger than physical nuclei. 
However, if semi-hard gluons are also included as sources 
of color charge (as they should be) then /i^ is defined as 

M 



1 / 1 

dxy—q{x,Q'^) 



m - 1 



gix,Q') (3.2) 



where q, g stand for the nucleon quark and gluon struc- 
ture functions at the resolution scale of Q of the physical 
process of interest. Using the HERA structure function 
data, Gyulassy and McLerran estimated that /i < 5 fm^^ 
for LHC energies and /i < 2.5 fm^^ at RHIC energies. 
Since semi-hard gluons appear at x values which corre- 
spond to rapidities greater than the central rapidity in 
nuclear collisions, we have to describe also the longitudi- 
nal dynamics of the fields in the collision. Such a descrip- 
tion is possible on a gauge lattice with a finite extension 
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into the longitudinal direction as used in the subsequent 
calculation. For a resolution of the longitudinal dynam- 
ics, the fact that the nuclei move with a speed smaller 
than speed of light becomes very important. This situa- 
tion is depicted in the following Fig. 1. 
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FIG. 1. The trajectories of two colliding nuclei (thin solid 
lines) are displayed together with the light cone (dot dashed 
lines). All components of the nuclei stay in the light cone 
throughout the collision since they propagate slower than 
speed of light. Consequently, the nuclei have a finite exten- 
sion in the z-direction. The dashed horizontal line indicates 
the direction in which we carry out the Fourier analysis of the 
initial state of a single nucleus. 

The finite longitudinal extension of the nuclei is also nec- 
essary for the propagation of the fields on the lattice. 
In the following, we demonstrate that the proposed de- 
scription allows to generate an "initial state" with a finite 
but small charge fluctuation. For the above chosen lat- 
tice constants each nucleon is covered by 40 lattice cells 
of the dual lattice. Each of these boxes is initially oc- 
cupied by an even number n^, of color charged particles 
such that the sum over all charges in a box is zero. The 
initial charge density thus is zero on a length scale of the 
lattice spacing and larger. At the initial time both 
nuclei, i.e. the particles, are boosted with the same 7 
but into opposite directions towards the center of colli- 
sion. The lattice constants a; and the time step width 
At are Lorentz transformed accordingly. The initial dis- 
tance do of the nuclei is chosen 16000 lattice spacings. do 
is Lorentz contracted to 45.07 fm in the center of veloc- 
ity frame. According to the above fixed extensions into 
transverse directions and size of the nuclei, we use a lat- 
tice that comprises a total of 4^ x 16384 points. We repeat 
the calculation for a longitudinally refined lattice which 
comprises a total of 4^ x 32768 points. In this case, the 



propagation starts at do = 32000 lattice spacings while 
the longitudinal lattice spacing is reduced to 03 = 0.1 fm, 
i.e. a single nucleus is then covered by 4 x 4 x 120 lattice 
cells. 

While the nuclei translate on the lattice, we adia- 
batically increase the coupling constant with a rate of 
Ag/At < (3000a) until g has reached its final value 
{g — 2.0). The propagating particles transport color 
charges from one box to another and the charge den- 
sity becomes different from zero and shows fluctuations 
in coordinate space. The spacial charge density on the 
lattice is defined by 



E 



(3.3) 



where (t) denotes the set of particle indices of the par- 
ticles which are in the dual box of the site x at time t. 
The upper index denotes the color while the lower index 
refers to the lattice site at which the density is evaluated. 
In Fig. 2, the distribution d(t^ q'^) of the charge per box 
is displayed at several different time steps t„. The charge 
with color c per box at site x is defined as 



(3.4) 



At time t = all q^{t — 0) are zero and d{t,q^) is 
different from zero only in the bin around q^ = 0. In 
the continuum limit a/ — > this would correspond to 
d{t = 0, g^) S{q^). From such a configuration we can 
generate a configuration possessing a finite charge fiuc- 
tuation or a finite /i respectively. 

Fig. 2 



tsoooo 



to 




-2 2 4 

charge per box [10"^] 



FIG. 2. The computed distributions of the charge per box 
for a boosted nucleus (7 = 106.5, 4 x 4 x 40 dual lattice cells, 
nt — 6, ai = a2 = as = 0.3 fm) for increasing times. The 
time is indicated at the curves by the time step numbers t„. 
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As already mentioned above, the initial positions, mo- 
menta and color charge orientations are chosen randomly. 
The chaotic nature of the Yang-Mills equations enhances 
strongly the convergence towards color distributions in 
position space which fluctuate in a highly statistical man- 
ner, i.e. the information about the initial distributions 
becomes strongly washed out after some time. 

Fig. 2 shows for one color, that the distribution of 
the charge density increases in width monotonously with 
time. The corresponding width is displayed in Fig. 5. 
Surprisingly, the shape of the distributions seems to be 
different from Gaussian distributions. This difference is 
explained in the following way. The shape of the dis- 
tributions in Fig. 2 results from a superposition of two 
Gaussian distributions of different widths which become 
equal after long times. In order to show that this ef- 
fect results from the Lorentz boost, we display in Fig. 
3 the color charge per box distribution of an unboosted 
nucleus. These distributions become Gaussian after rel- 
atively short times as indicated by the small time step 
numbers i„ in the figure. 



the direction of the boundary transition. On a Lorentz 
contracted lattice, charges that cross transverse bound- 
aries experience a much stronger rotation than charges 
which cross longitudinal boundaries because the trans- 
verse links are longer and the transverse color electric 
fields are stronger than the longitudinal components. In- 
deed, for small times we find a shape of the charge dis- 
tribution similar to that in Fig. 2 in the time evolution 
of an unboosted collection of particles on a (by arbitrary 
choice) deformed lattice. 

The curves in Fig. 2 from iioooo to t4oooo contain a 
slightly broadened high "peak" in the center and a wide 
part at the bottom. In order to show that the slow 
increase of the width of the central peak results from 
longitudinal color rotations we refine the lattice in the 
longidutinal direction by a factor three as already men- 
tioned above while keeping the volume of the nucleus 
constant. The number of particles increases by a factor 
three while the length of the charge vector is reduced by 
a factor three. 



Fig. 3 
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FIG. 3. The same distributions as in Fig. 2 but for an 
unboosted nucleus (10 x 10 x 10 dual lattice cells, nt — 6, 
a\ = a2 = a3 — 0.3 fm). 

For a boosted nucleus, the lattice is Lorentz contracted 
in the longitudinal direction. The randomly generated 
coordinates of the particles imply a color charge fluctua- 
tion which appears after short times in the distributions. 
In addition to the propagation of the charges, a rota- 
tion in color space occurs at a charge when it crosses 
a cell boundary. These rotations lead to an additional 
smoother broadening of the charge distributions in each 
color. The rotation depends on the length of the link 
and the strength of the color electric field pointing in 



Fig. 4 
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FIG. 4. The same distributions as in Fig. 2 but for a refined 
lattice (7 — 106.5, 4 x 4 x 120 dual lattice cells, nt — 6, 
ai = a2 = 0.3 fm, 03 — 0.1 fm). The dashed curve shows the 
Gaussian fit of the distribution at time step tsoooo- 

In Fig. 4, we display the color charge distributions for 
a boosted nucleus in the case of a finer longitudinal lat- 
tice. The central peak decays gradually in time but stays 
narrow in comparison with the wide distribution at the 
bottom and in comparison with the central peak in Fig. 
2. The distribution displayed for the time step ^5000 re- 
sults essentially from charge propagation. The width of 
this distribution corresponds exactly to the length qo of 
a charge vector. At later times color rotation seems to 
smear out the distribution. At ^50000 which is close before 
the collision, the distribution has a Gaussian shape. In 
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principle one would have to evolve much longer in time 
until the central peak is smeared out. The time step 
numbers in Fig. 4 are much larger (by a factor 7^) than 
in Fig. 3. This is explained by the fact that the boosted 
particles propagate essentially into the longitudinal di- 
rection. In addition, the time step width is reduced by a 
factor 7~^. 

In the subsequent calculations, we neglect the resulting 
final changes of the width because they are small. We find 
the same time behavior of the distributions for all colors. 
The fluctuation of the charge density (JpLatt associated 
with the width of the distributions shown in the figures 
Fig. 2 and Fig. 4. Numerically, (JpLatt is determined by 

ispLuf--=j^uj:^^r^ (3.5) 

where we make use of the fact that 

EE'?^' = 0- (3.6) 



Fig. 5 
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FIG. 5. The fluctuation of the charge density 5pLatt on 
the lattice (upper solid curve) is displayed for a Lorentz con- 
tracted nucleus (see units on axis of ordinates) as a function 
of time over the number of time steps. The time step width 
is At — 0.03/7 fm. The dashed curve shows the same but 
averaged over two points. 

In Fig. 5, we display (JpLatt ^f the distributions shown 
in Fig. 2 as a function of time. Fig. 5 shows, how 
^P'iatt evolves after the Lorentz boost when the nucleus 
propagates along the contracted lattice. The curves show 
that ^p£att begins to converge at large times. In Fig. 6, 
we display ^p£att '^^ ^^e distributions shown in Fig. 4 as 
a function of time. The curves are plotted for all colors. 



We now make an estimate of the transverse average 
charge density /i used in the McLerran Venugopalan 
model jisf . The correct way to determine this quantity 
is rather involved since /i is related to a surface charge 
density in the transverse plane. The calculation of this 
surface charge density requires a parallel transport on 
the SU(2) mainfold of all charges in the longitudinal di- 
rection from the spacial site of the charge into a fixed 
transverse plane, e.g. the central transverse plane of a 
nucleus. The surface charge density which follows from 
that procedure leads to a distribution of charges per pla- 
quette in the transverse plane. For reasonable statistics, 
the transverse extension of the lattice has to be large 
enough. The transverse extension of our lattice how- 
ever is small and we therefore give a simple estimate. 
We assume that, according to translational symmetry of 
the unboosted three dimensional lattice, the distributions 
shown in Fig. 2 and Fig. 3 are equal to the distribution 
which we would obtain for each transverse layer of boxes 
from a lattice with large transverse extension. The con- 
volution over all transverse charge density distributions 
leads to the surface charge density distribution in one 
transverse plane. Supposing that the distributions have 
a Gaussian shape, we determine /i by 

M = ^2^_ I E ^P'iattaj.\faia^3. (3.7) 

^ c 

iVa denotes the number of lattice cells which cover a nu- 
cleus in the longitudinal direction and corresponds to the 
number of convolutions. 



Fig. 6 
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FIG. 6. Same as Fig. 5 but for a lattice with 4 x 4 x 2^^^ 
points and constants ai = a2 = 0.3 fm, as = 0.1 fm. A nucleus 
is covered by 4 x 4 x 120 dual lattice cells. 

In the calculation of the curves in Fig. 5 and Fig. 6, we 
have used = 6. According to Fig. 5 and Fig. 6, we 
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obtain the final values /i = 0.06 fm~^ and fj, = 0.04 fm~^. 
These values stay below the estimate of Gyulassy and 
McLerran |Q. The value of fi or 6p respectively is 
one important criterion for the construction of the initial 
state before the collision. Since the charge Qo — ^3/4 is 
fixed for SU(2) quarks and g — 2 according to [Q, rib and 
the transverse lattice constants ai, 02 are the remaining 
parameters which determine fj,. 

Another important criterion is that the distribution of 
soft gluons and semi-hard gluons should reproduce the 
longitudinal momentum distribution g{x, Q^) which is 
observed experimentally. To reproduce g{x, Q^) would be 
far too ambitios but at least we have to consider the mo- 
mentum distribution. The gluon distributions in momen- 
tum space are simply related to the Fourier transformed 
gauge field A'^{k) by \2Ti {Au{k)Au{k))\ (no summation 
over v. The light cone model of McLerran and collabora- 
tors describes the initial state through pure gauge fields 



= 0, 

A' = e(z-)Q{-z+)a\{x,y) + 

eiz+)e{-z-)al{x,y) 1 = 1,2. 



(3.8) 
(3.9) 



The coordinates z^ are defined as z^ = ^{t ± z). 
In the following, we call the direction I = 3 which is 
parallel to the collision axis the "longitudinal" direction 
and E, A-field amplitudes associated with longitudinal 
links are called the "longitudinal fields" i?L- The fields 
Et in transverse directions are called "transverse" . The 
Fourier transformation of the field of one nucleus along 
the dashed line in Fig. 1 leads to a distribution of the 
form |Ai(fc3)P ~ 1/^3- The electric field follows from the 
equations ( |3.8| ) and (^) by the definition E = ~dtA. 
The transverse electric field along the dashed line Fig. 
1 has the z-dependence Ei ~ 5{z — zq) {i = 1, 2) while 
the longitudinal electric field Ez = ot E^ = respec- 
tively. The Fourier spectrum of the color electric field is 
therefore a constant and shows no fca-dependence. The 
following results presented in this section have been cal- 
culated for a single nucleus which propagates on a lattice 
of the size 4^ x 16384 points. In Fig. 7, we display 
as a function of the Fourier index k^, = k'ia^/{2'K) at the 
time step isoooo j the longitudinal Fourier spectrum of the 
transverse gauge field which we calculate according to 



- 2n,n,In?-i) ^ 5: E ^?(^' ^3) 



x,y c /=1,2 



(3.10) 



The corresponding Fourier spectrum of the transverse 
color electric field is calculated as 



/f^(i,fc3) 



(3.11) 



at the same time step. The gauge fields Af(i, x) have 
been determined on the lattice through time integration 



of the color electric fields along the whole time evolu- 
tion. The momentum distribution in Fig. 7 agrees (if 
one compares the normalized distributions) with the re- 

suit obtained for fj, {t^k^). This check tells us that 
the gauge fields Ai{t,x) displayed in Fig. 9 (a) are true 
in a sense that they really correspond to the color elec- 
tric fields. In Fig. 8, the Fourier spectrum of the lon- 
gitudinal gauge field is displayed on a logarithmic scale, 
i.e. \ogiQ{\ f^\t,k^)\'^k1) is plotted as a function of the 
Fourier index k^ = k^a^^ / {27:) . The figure shows that 
the longitudinal components of the fields carry very low 
longitudinal momenta in comparison to the longitudinal 
momenta in the transverse components. Most of the lon- 
gitudinal field energy is in modes with fc < 34 while the 
transverse field energy is in modes up to fc = 3400. 
Fig. 9 (a) shows for all colors the gauge fields generated 
by a nucleus which moves from the left to the right in the 
figure. The steep front of the field distribution indicates 
the position of the nucleus. 

Fig. 7 




FIG. 7. Longitudinal Fourier-spectrum of the transverse 
gauge fields [/^^ (isoooo, '^s)^^! in units of 10^ fm"'' plotted 
as a function of the Fourier index kz. 

Lets assume that the color electric field Ef{z) is finite 
and constant inside of the nucleus and zero outside of the 
nucleus, i.e. Ef{z) = 1 for ^ < z < and Ef{z) = 
for z > ^ or z < The Fourier transformation of 

such a static rectangular field distribution has the form 



23/2 sin(^fc3) 



(3.12) 



from which we conclude that the half period of \Ef{k^)\'^ 
in fc-space is Afcs = This difference can also be 

expressed in terms of the Fourier index as Afc = . 
For a moving nucleus, the Fourier spectrum is shifted as 
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can be seen in Fig. 7. The width of the large hump in 
Fig. 7 corresponds to the longitudinal extension of the 
boosted nucleus which is Az = 40a3 while N3 = 16384. 

Fig. 8 




FIG. 8. The longitudinal Fourier-spectrum of the longitu- 
dinal gauge fields is shown as a function of the Fourier index 
ks on the logarithimic scale as 21ogjQ (^1/^^' (fgOOOO, fca) Ifcs) . 
The argument of the logarithm is given in units of fm~^. 

In Fig. 9 (a), we display a snap shot of the gauge field 
amplitudes \A {t,z)\'^ taken as average over the trans- 
verse planes at the time step iiooooo- Fig. 9 (b) shows 
for the same time step the corresponding transverse en- 
ergy density distribution which we define as 



{z) = Tr{£i£i)+Tr{£2£2)- 



(3.13) 



Further below, we will also need the longitudinal energy 
density distribution which is defined 



(3.14) 



The figures Fig. 9 (a) and Fig. 9 (b) clearly show that 
the gauge field behind the moving color sources is mostly 
a pure gauge field because the amplitude of the color 
electric field behind the nucleus is zero. Large amplitudes 
of the color electric field appear only inside of a nucleus 
which has a finite extension in the longitudinal direction. 
The gauge fields displayed in Fig. 9 (a) correspond to 
the gauge fields of the light cone model according to the 
ansatz in Eq. (3.S). In the light cone model, behind 
a nucleus, they are constant along lines parallel to the 
collision axis. These fields, however, are not unique due 
to gauge freedom. Therefore, the gauge fields in Fig. 9 
(a) can be very different from the ansatz in Eq. ( p.§| ). 




Fig. 9 



FIG. 9. The transverse gauge fields \A (iiooooojz)^ aver- 
aged over the transverse plane and squared are displayed for 
each color in the upper panel. In the lower panel, we show 
the corresponding transverse color electric field energy density 
Wt (iiooooo,^)- 



IV. THE COLLISION 

In the previous section, we have outlined how a collision 
can be simulated on a gauge lattice with small trans- 
verse extension and how the initial state can be gen- 
erated within the model. In this section, we focus on 
the time period after the collision and discuss results ob- 
tained from simulations of a collision on two different 
lattices. In one case we performed a calculation on a lat- 
tice of the size 4x4x2^^ points using lattice constants 
fli = 02 = 0.3 fm and = 0.1 fm. Accordingly, a single 
nucleus is covered by 4 x 4 x 120 dual lattice cells. In 
the following the refer to this paramerization as the "case 

(1) " . In the second case, we performed a calculation on a 
lattice of the size 4 x 4 x 2^^* points using lattice constants 
ai — 0.3 fm(/ — 1,2,3). A single nucleus is thus covered 
by 4 X 4 X 40 lattice cells. We refer to this parametrization 
as the "case (2)". A time-step width of At = 0.03/7fm 
is used in both cases. In the case (2), we construct initial 
states as "realistic" as possible within the model. We 
try to ajust the glue field energy of the nucleon, to ob- 
tain a color charge fluctuation fi in the range estimated 
by other authors and to describe the longitudinal size of 
the Lorentz boosted nucleus. We expect that this calcu- 
lation provides a non-perturbative estimate for soft glue 
field scattering and glue field radiation at high energies. 
Due to the high color field energies the simulation of case 

(2) is close to the limit of acceptable numerical precision. 
An improvement requires finer but also larger lattices ex- 
ceeding our computational resources. Instead, we refine 



9 



the lattice in case (1) keeping the size small. In this case 
the fields carry less energy before the collision but the res- 
olution of the field dynamics is higher and more precise. 
Case (1) therefore provides a control over the qualitative 
nature of the results while we expect some quantitative 
results from case (2). Most of the figures presented below 
are shown for case (1). The curves agree in a qualitative 
manner with results obtained in case (2) but they show 
less intensity. 

In Fig. 9 (b) of the previous section, we have shown 
the transverse color electric field energy density distribu- 
tion (z) of a propagating nucleus for one time step. 
To get a full picture of the evolution of {t, z) and 

{t, z) during the generation of the initial state when 
the two nuclei approach each other, we display in Fig. 
10 (a) and (b) the energy densities of the transverse and 
longitudinal components of the i?- field every 10000 time 
steps for the case (1). 
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FIG. 10. The transverse (a) and longitudinal (b) energy 
densities of the color electric fields are displayed for two ap- 
proaching nuclei. The numbers on top of the distributions 
indicate the number of steps in the time integration (on top 
of left nucleus) or time in units of fm (on top of right nucleus) . 

We verified, that the corresponding B-field energy densi- 
ties w'^^^ (z) are almost identical with w'r^^ [z). After 1000 
time steps (not shown in the figure) the transverse _E-field 

i E^ 

energy Wr^ is by about a factor 7/2 larger than the lon- 

( F1 

gitudinal i?- field energy . This ratio increases to 7 
at 50000 time steps and remains constant afterwards. 
As Fig. 10 (a),(b) shows, two configurations representing 
two nuclei propagate remarkably stably (nucleus 1 from 
left to right and nucleus 2 from right to left) over long 



distances on the lattice. The collision begins at ^53333 
where the fields of both nuclei make first contact. 
In Fig. 11 at time step tsoooo of case (1), the particle 
density (in units of 20 on the ordinate), which is defined 
as the number of particles per bin of width a/10 in longi- 
tudinal direction, is superimposed on the transverse and 
longitudinal field energy densities (multiplied by a factor 
10 in the figure) for comparison. The particle density 
defines the extension of the nucleus in the longitudinal 
direction. The mismatch between particle density and 
field energy density decreases when the density of lattice 
points covering the nucleus in the longitudinal direction 
is increased. This requires a shift of the cutoff (here at 
kc — 2.0 GeV for a — 0.3 fm) to higher momenta. 

Fig. 11 
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FIG. 11. Comparison of the transverse (T) and longitu- 
dinal (L) color electric field energy distribution in a single 
nucleus with the distribution of the particles. The longitu- 
dinal energy density has been multiplied by a factor 10 for 
better comparison. The particle density (P) has been divided 
by a factor 20. 

The tails of the field energy densities result not only 
from lattice dispersion effects. The lattice is Lorentz con- 
tracted in the longitudinal direction and particles have 
large longitudinal momentum components. Therefore, 
particles primarily transfer energy into longitudinal links. 
The nonlinearity of the Yang Mills equations provides 
a mechanism transferring energy from longitudinal into 
transverse degrees of freedom. Field amplitudes which 
are left behind on the longitudinal links are canceled by 
following particles of rotated color charge. 



4 

Fig. 10 



10 




Fig. 12 



z— Coordinate [fm] 

FIG. 12. The transverse color electric field energy density 
is displayed for case (1) at various different time steps after 
the begin of the collision. The total time step number t„ and 
the relative time t counted from the first contact of the nuclei 
is indicated in each panel. The collision begins at the total 
time step fsssss or at the relative time t = Ofm respectively. 
The particle distributions (multiplied by a factor 1/10) are 
included in panel (a) and (b). 

Fig. 12 displays snap shots of {tn,z) 1667, 6667, 
11667, 16667, 21667 and 26667 time steps after the be- 
gin of the collision. In the panel (a) and (b) the particle 
distributions in coordinate space (multiplied by a factor 
1/10) are displayed together with the field energy distri- 
butions to indicate the position of the receding nuclei. 
As compared to Fig. 11, the width of the distributions 
Wrp (t, z) moving with the particles is considerably in- 
creased (about 30%). This sudden change of the width 
occurs during the overlap of the nuclei. At the same time 
the field amplitudes at the position of the particles have 
decreased. We would expect that this leads to an in- 
creased energy transfer from the particles into the fields. 
To explore this behavior, we define the electric field en- 
ergy for each color 



aia2a3 



xeX 1 = 1 



(4.1) 



Note, that we do not sum over the color index c. The 
magnetic field energy Wr^ (t) is defined in the same man- 
ner. In the calculations, we find tha.tWr'it) ~ wr'{t). 

We calculate the integrated field energy as a 

function of time. Fig. 12 will be further discussed below. 
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FIG. 13. Same as in Fig. 12 but for the longitudinal color 
electric field energy density w'^\z). 

{ E) 

In Fig. 14, we display for case (2) Wc (t) for each 
c in a time interval around the time of first contact 
isoooo- In the time region between isoooo and igoooo, 
we observe a kink (it is a bit washed out) in Wc'^\t), 
which increases remarkably in the collision (see also Fig. 
15). In the upper left corner, we display for case (1) 
W^^\t) — J2c Wc^\t) in a time interval around ^53333. 
The angle between the dashed line and the dot-dashed 
line indicates the kink. The kink in Fig. 14 corresponds 
to an increased energy transfer from the particles into 
the fields while the nuclei recede. This tells us that the 
energy deposit between the receding nuclei results not 
only from the interaction between the fields during the 
overlap time but glue field radiation contributes. The 
contributions to glue field radiation appear in the longi- 
tudinal field energy density distributions shown in Fig. 
13. Energy is transfered from the particles into the fields 
through longitudinal links. Before the collision during 
the generation of the initial state this transfer occurs 
smoothly generating coherent fields. The strong increase 

( E\ 

of Wj^ {t, z), at times larger than the time of the kink and 
further away from the center of collision, results from the 
fact that more energy is transfered into longitudinal links 
than the closest transverse links can absorb from them 
during the short time in which a nucleus passes a cer- 

( E) 

tain longitudinal link. The distribution wf^ (z) increases 
sharply for collision times larger 0.3 fm and arrives at a 
first maximum at i ~ 0.5 fm. The next minima appear 
at a time t ~ 0.7 fm which could presumably be inter- 
preted as a thermalization time for the soft glue fields. It 
corresponds to the FWHM of the valley in Fig. 13 (c). 



11 



Fig. 14 




7 8 9 

Number of Time Steps in Units of [1 O ] 



FIG. 14. The total color electric field energy W^i^'(t) for 
each color is displayed for case (2) in the time interval between 
^60000 and tiooooo • The panel in the upper left corner displays 
the total electric field enery for case (1). The lines indicate a 
kink. 

When the two nuclei pass through each other, the parti- 
cles experience the combined fields of both nuclei. The 
corresponding chan ges of the field amplitudes enter into 
the r.h.s. of (2.3) and modify the orientation of the 
color charge vectors gf(t) during the rather short over- 
lap time of 0.11 fm. This induces net color charge cur- 
rents resulting in radiation of gluons. Since the longitu- 
dinal momenta of the particles are large compared with 
their transverse components, color charges move essen- 
tially parallel to longitudinal links and induce electric 
fields on these links. Before the collision, the charge of 
the following particles was polarized such that these fields 
were canceled. After the collision this is no longer true 
over a long period of time. A continuation of the time 
evolution has shown that the longitudinal energy density 
of the fields remains as large as in Fig. 13 and decreases 
slowly at very large times. 

The excess of energy remaining in the longitudinal links 
propagates then into transverse directions. These longi- 
tudinal fields possess only transverse momenta and re- 
sult in gluon radiation into transverse directions. This 
mechanism describes gluon radiation from the particles. 
The fact that the transverse extension of our lattice is 
very small and periodically closed does not really allow 
an expansion of the radiated glue fields into transverse 
directions. They collide with other radiated gluons and 
are stopped. This feature is somewhat unphysical and 
can only be removed by taking lattices with large trans- 
verse extensions. The radiated glue fields remain there- 
fore essentially on the longitudinal links and decay slowly. 
However, at least it allows one to determine how much is 



radiated from a certain volume. 

We therefore define the color separated longitudinal elec- 
tric field energy 



aia2Q3 \ - 
2 ^ 



(4.2) 



which is a function of time. Fig. 15 displays W£ c (^) ^'^^ 
all colors. The kink in the time dependent field energy 
of Fig. 14 appears much stronger in Fig. 15. The curves 
in Fig. 15 show essentially the cumulative time integral 
of how much energy is radiated from the particles at in- 
creasing time. 



Fig. 15 
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FIG. 15. Longitudinal color electric field energy as function 
of time. 

In the following, we discuss the contribution to the energy 
deposit which results from the interaction of the fields. 
The contribution of the particles will be discussed later. 
We discuss the results for two choices of the lattice size 
and lattice spacing. The numbers, obtained in the case 
where we use a 4 x 4 x 16768 lattice are given in paren- 
theses. As we see from Fig. 12, a considerable fraction of 
the transverse field energy is deposited around the center 
of collision between the receding nuclei. 

At the time step iieoooo in the case (2), about 440 GeV 
of the total transverse electric field energy Wr^^"* = 
2810 GeV are left in the spatial region between zi — 
3.08 fm and Z2 = 43.08 fm (the center of collision is at 
z = 23.08 fm). The field energy density between the nu- 
clei is about a factor 1/100 smaller compared to the field 
energy density at the position of the particles in the re- 
ceding nuclei. Due to the large space volume between 
zi and Z2 however, the energy deposit amounts to about 
16% of the total energy in soft field modes. 
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Fig. 13 shows that the longitudinal energy densities 
are small but finite around the center of collision. The 
calculation shows in case (2) that within short times 
At — 1.0 fm after the begin of the collision about 23% 
(7.42 GeV) of the total field energy deposit (31.98 GcV) 
that is left in the region between zi — 22.08 fm and 
Z2 = 24.08 fm are carried by longitudinal components of 
the £'-field. The average ratio /wi^^ in this region 
is about 0.6 which is large compared to the above men- 
tioned factor 1/7. This ratio appears to be even larger 
in case (1). Within the short time At — 0.3 fm after 
the begin of the collision about 29% (664 MeV) of the 
total field energy deposit (2270 MeV) left in the region 
between zi = 15.08 fm and Z2 — 15.68 fm are contained 
in longitudinal component of the i?-field. 

Since the color charge and color current density is zero, 
the linearized Yang-Mills equations (2.4) are homoge- 
neous in this region and allow only for solutions with non- 
zero amplitudes into transverse directions in relation to 
the direction of their energy flow. Figure 10 displays the 
fraction with momenta pointing into transverse direction 
relative to the coUison axis. The longitudinal field energy 
density in the narrow space region between zi = 15.08 fm 
and Z2 = 15.68 fm is not a result of gluon radiation. A 
calculation with the right-hand side of (2.2) and ( |2.3| ) set 
to zero just before the nuclei make first contact yields 
the same results in that narrow region, which finds its 
explanation in a transfer of energy from the propagat- 
ing transverse fields to the longitudinal fiel ds d ue to the 
non-linear coupling between £^t and in (2.4). Also in 
collisions of Yang-Mills wave packets, where no particles 
are taken into account, we observed a strong increase of 
the longitudinal fields around the center of collision dur- 
ing the overlap of the wave packets |l^ . 
When the two nuclei overlap during the collision, the 
fiel ds a re superposed. As a result of the nonlinear terms 
in (2.4), which act as source terms for the longitudinal 
fields, the amplitudes start to grow rapidly during the 
short overlap time (of the particle clouds) of 0.11 fm. This 
mechanism describes the scattering of soft gluons which 
occurs over a much longer period of time. To analyze this 
process within the model, we define the following quan- 
tities. The Foynting vector in the adjoint representation 
is denoted as 



S:=c£x B. 



(4.3) 



With t he t ransverse and longitudinal components of the 
vector (||^) we define the total transverse and longitudi- 
nal energy currents 



^70000- Fig. 16 (a) shows that shortly after the first con- 
tact of the nuclei there is a sudden increase of the trans- 
verse energy current in the fields. This increase describes 
soft gluon scattering into transverse directions. The slope 
of the dot-dashed line indicates the current correspond- 
ing to the energy rate transfered from the particles to the 
fields if no collision would occur. The dot-dashed line rp- 
resents an underground which has to be subtracted from 
the curve to determine the contribution that results from 
scattering. A very similar behavior has recently been 
found in collisions of color polarized wave packets on the 
three-dimensional lattice ||l^. At the same time, Fig. 16 
(b) shows that the longitudinal energy current decreases 
remarkably during the overlap of the fields. The width of 
the inverse peak in Fig. 16 (b) corresponds to a time of 
0.18 fm which is the overlap time of the fields. After the 
overlap, the curve approaches the dot-dashed line almost 
as close as before the collision. This leads to the conclu- 
sion that the inverse peak is a interference phenomenon. 
Nevertheless, it determines the time interval in which 
the fields interact. The steady increase of the longitudi- 
nal current before the collision results from the energy 
transfer from particles into fields. In a similar study for 
case (2), we find that the transverse energy current keeps 
increasing much stronger long after the collision. These 
additional contributions describe glue field radiation. 
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In the following Fig. 16, we display irit) (panel (a)) and 
ihit) (panel (b)) for the time interval between iaoooo and 



FIG. 16. The transverse (panel (a)) and longitudinal (panel 
(b)) energy currents irit) and iL{t) are displayed in a time 
interval around the beginning of the collision. 

To give an estimate for glue field radiation, wc inte- 
grate the curve in Fig. 13 (f) from zi = 8.38 fm to 
Z2 = 22.38 fm and find an energy deposit of 64.62 GeV. 
This energy divided by the volume covered by both 
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particle distributions and divided by the time oi t = 
7.51 fm since first contact yields a radiation power of 
about 115GeV/fm*. In the case (2) we find the re- 
sult llOGeV/fm'* over a time of 23 fm. We now give 
estimates for contributions from soft scattering and ra- 
diation. The total field energy at the time step tsoooo 
amounts to 757.24 GeV in case (1) while the field en- 
ergy deposit in the space region between zi = 8.38 fm 
and Z2 = 22.38 fm is 94.74 fm. We conclude that 8% of 
the total field energy are deposited between the nuclei at 
t = 7.51 fm. In case (2), we find 26% at t = 23 fm. 
Finally, we analyze the momentum distribution in the 
transverse color electric fields. In the previous section we 
have discussed the Fourier spectrum of the initial state 
of a single nucleus. Now, we transform the fields of two 
nuclei propagating at the same time on a lattice of the 
size 4 X 4 X 2^*^ according to case (1). In the upper panel 
(a) of Fig. 17, we display the normalized longitudinal 
momentum distribution of two nuclei at time step isoooo 
shortly before the beginning of the collision. With Eq. 
(3.11) it is defined as 



FIG. 17. The normalizd Fourier spectra of the transverse 
color electric field before (a) and after the collision (b). The 
snap shots are taken at the time steps tn indicated on top of 
the curves. In the upper left corner of panel (b)), we display 
the cumulative integrals over the spectra. 

In the upper right corner of Fig. 17 (b), we display the 
cumulative integrals over the distributions 



(4.7) 



fT\t,ks) 



dk3f^^\t,k3) 



(4.6) 



In the lower panel (b), we display the momentum distri- 
bution after the collision at the time step tsoooo- A com- 
parison of Fig. 17 (a) with Fig. 17 (b) shows evidence 
for the production of a final state with a momentum dis- 
tribution that dramatically differs from that of the initial 
state. 



Fig. 




We compare the results for isoooo (solid curve) and for 
^80000 (dashed curve). Above k = 2200, the dashed 
curve is clearly below the solid curve which shows that 
energy is transfered from high frequent modes into low 
frequent modes during the collision. This behavior has 
also been found in collisions of color polarized Yang-Mills 
wave packets in 1-1-1 dimensions ^ and later in collisions 
of Yang-Mills wave packets in 3-1-1 dimensions JM| and 
could presumably describe particle production ||32| |. Our 
calculation shows that this pure classical phenomenon oc- 
curs also in collisions of Yang-Mills fields which are not 
polarized in color space. 



V. SUMMARY AND CONCLUSION 

In the framework of a classical effective approach of QCD, 
color charged particles have been employed to generate 
the color source effects of the valence quarks in ultra- 
relativistic nuclear collisions. Classical Yang-Mills fields 
.,7 have been used to describe soft modes of the gluons. 
We have studied the time-evolution of such a system 
of particles and fields by solving the coupled system of 
Wong equations and Yang-Mills equations, as a transport 
model for the reaction of nuclei in collisions at high en- 
ergies. We have carried out numerical simulations of the 
collision on a long three dimensional gauge lattice with 
small transverse extensions. Results have been presented 
which demonstrate the properties of such a description 
from various aspects and provide estimates for contribu- 
tions at small transverse momenta. We have discussed 
the similarities and the basic differences of the approach 
presented in this paper as compared to the light cone 
source model |ljjl^^. It has been demonstrated that 
a "initial state" can be generated which describes the lon- 
gitudinal extension, the color charge fluctuation and the 
glue field energy of the nuclei before the collision. The 
shape of the longitudinal momentum distribution of the 
soft gluons in the initial state as obtained in our calcu- 
lation is similar to the shape of the experimentally know 
gluon distribution at small x. 

As a function of time, we have calculated the transverse 
and longitudinal energy densities of the color electric 
fields. For an effective coupling oi g = 2 and for an energy 
of lOOGeV/u, our calculation predicts that a relatively 
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large fraction (over 8%) of the total soft field energy (at 
small x) contributes sX t = 7.5 fm after the first contact 
to the energy deposit between ± 7 fm around the center 
of collision zq- 

The colliding color fields cause soft gluon scattering. This 
has been demonstrated in a calculation of the transverse 
and the longitudinal energy current as a function of time 
throughout the collision. 

The non-linear nature of the Yang-Mills equations causes 
a change of the color charge distribution in the two col- 
liding clouds of particles during the overlap time which 
results in an increase of the charge fluctuation. The re- 
sulting color charge currents lead to an averaged gluon ra- 
diation from the particles which amounts to 110 GeV/fm'' 
(radiation power per particle volume) over times of 7 fm 
and even 20 fm. 

The longitudinal momentum distributions of the trans- 
verse color electric fields have been calculated. The fields 
were not polarized in color space. Our results show that 
an energy transfer from high frequenzy modes into low 
frequenzy modes occurs in the collision. The same obser- 
vation has been made by Hu et al. Q studying collisions 
of color polarized pure Yang-Mills wave packets on a one- 
dimensional SU(2) gauge lattice. The Fourier analysis of 
section 5 shows that the existence of this phenomenon 
does not require color polarization in the initial condi- 
tions, e.g. before the collision. 

In this paper we have discussed a new approach to de- 
scribe ultra-relativistic heavy-ion collisions in a coher- 
ent way in 3-1-1 dimensions. We have presented results 
from numerical model simulations to explore the possi- 
bilities for describing results from experiments. A com- 
bination of our approach with parton cascades and ultra- 
relativistic quantum molecular dynamics could allow to 
trace further back from experimental data into the early 
stage of the collision. The approach for itself could de- 
scribe the formation of the quark gluon plasma within 
the classical mean field approximation. 
At the present stage, we are still far from accomplishing 
such goals. It would be important to increase the ex- 
tension of the lattice into transverse directions to cover 
two nuclei completely. This would allow to determine 
the soft transverse momentum distributions and the time 
evolution of single transverse modes. Also a study of the 
dynamics of the particles will be necessary. This requires 
that the particles are fully coupled to the fields. The mo- 
mentum distributions and rapidity distributions of the 
particles could be compared with experiment. 
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